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Abstract 

A closed form expression for the ground state energy density of the general 
extensive many-body problem is given in terms of the Lanczos tri-diagonal 
form of the Hamiltonian. Given the general expressions of the diagonal and 
off-diagonal elements of the Hamiltonian Lanczos matrix, a„(A^) and Pni^), 
asymptotic forms a{z) and (3{z) can be defined in terms of a new parameter 
z = n/N (n is the Lanczos iteration and N is the size of the system). By 
application of theorems on the zeros of orthogonal polynomials we find the 
ground-state energy density in the bulk limit to be given in general by £q = 
inf [a{z) -2p{z)]. 
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I. INTRODUCTION 



Finding analytic means of calculating details of the energy spectrum of strongly inter- 
acting many-body systems has long been a goal of theoretical physics, with applications 
from low energy scales in condensed matter physics to high energy particle physics. In 
particular, accurate calculations of the ground state energy are important to settle issues 
relating to the character of the ground state, or vacuum state, itself. While there has arisen 
a range of theoretical tools to tackle problems that cannot be treated by perturbation theory, 
many suffer from the effects of uncontrolled approximation to the extent that each, alone 
cannot be trusted to give a reliable answer. It is only when several independent methods 
point in the same direction, that one can confidently state something about the character 
of a system. Examples include mean field theories neglecting correlations in fluctuations 
thereby reinforcing a tendency to order; exact diagonalisation and Monte-Carlo studies, at 
zero temperature, on relatively small clusters and the extrapolation to the bulk limit; fi- 
nite temperature Monte-Carlo simulation of somewhat larger systems, but at temperatures 
above the interesting low energy scales. 

In this paper a simple fundamental relationship is found between the tri-diagonal form 
of the Hamiltonian and the ground state energy density. A general theorem is proved for the 
ground state energy density in terms of the coefficients generated by the Lanczos method, 
evaluated in a limiting process incorporating convergence of the Lanczos iterates and the 
thermodynamic limit. The Lanczos method is based on the following recursion: starting 
from an appropriate trial state, one recursively generates new basis states with repeated 
application of the Hamiltonian 



Pn-l L J 

where q;„ = {vn\H\vn) and /3„ = {vn+i\H\vn) ■ At the nth iteration of this recursion the 
Hamiltonian matrix in the new basis of states is given by the tri-diagonal, T^, i.e. 
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The power of the Lanczos phenomenon is that the dimension of the tri-diagonal basis, 
required to describe the low-lying states of the system, determined by the recursion level 
is significantly smaller than the original basis. The outermost eigenvalues of T„ rapidly 
converge to those of the Hamiltonian. In computing the Lanczos coefficients, and 
exactly one is usually restricted to an early termination of the recursion in the analytic case 
treating large systems, or to small systems in a numerical calculation taken to complete 
convergence. To date all methods, however the Lanczos basis has been generated, require 
the numerical diagonalisation of the tri-diagonal T„ matrices for the ground state energy. In 
this paper we will demonstrate how the diagonalisation can be carried out analytically for 
extensive systems, thereby providing a solution for the ground state energy density in terms 
of the tri-diagonal form. By introducing a new parameter, z = n/N, where is the size of 
the system, the asymptotic forms a{z) and (3{z) are related to the ground state energy by 



II. THERMODYNAMIC LIMIT OF THE TRI-DIAGONAL FORM 

Using the initial state, \vi), one forms Hamiltonian moments, and from these the con- 
nected moments. 



The connected moments encapsulate the essential physics of the system because they scale 
with the size of the system, N, as 



So = inf [a{z) - 2 p{z)] . 



(3) 




(4) 
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(i/")e = C„iV. 



(5) 



Here is a quantitative measure of the size of the system whether it be the number of 
sites in a lattice of locahsed spins, or the volume in a continuum model with itinerant 
particles. Although this form is restricted to the ground state or vacuum sector of the 
model, generalisations can be easily made to excited states. 

Al/N expansion of the Lanczos matrix elements, an{N) and /3„(A^), reveals a surprising 
analytic property - a simple polynomial n dependence. In terms of the connected coefficients 
c„, this cluster expansion for the Hamiltonian density are |Jl|,0| 



ttn(iV) = ci + (n - 1) 
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We first define a parameter z = n/N which remains finite as the Lanczos recursion 
proceeds in the bulk limit, n — 00 and N 00. In this limit the expansions become a 



series m z, i.e. 



a{z) = lim an{N) = Ci + z 
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The cluster expansion guides us to the observation that, more generally, we have for 
the exact problem the asymptotic forms in the — > 00 and N ^ 00 regime 



an{N) = a{z) + 0{l/N), 
f3l{N)=(3'{z)+0{l/N). 



(10) 



We observe here a confluence of the two hmiting regimes: that of convergence of the Lanczos 
iterates and of the thermodynamic hmit into a single scaled Lanczos iteration number z. 

III. ORTHOGONAL POLYNOMIALS AND VAN DOORN'S THEOREM 

The connection between the Lanczos tri-diagonal form of the Hamiltonian and the 
associated system of orthogonal polynomials is simple. The characteristic polynomial 
Dn{x) = det(T„ — xin) of the Lanczos tri-diagonal matrix representation satisfies the fol- 
lowing recursion relation: 



which in turn defines Pn{x) = (— 1)"-D„(x) as an orthogonal polynomial, the zeros of which 
are the eigenvalues of the T„ matrices. The analytic forms for the an{N) and f3n{N) define 
a special class of orthogonal polynomials relevant to the many-body problem which are 
distinguished by a certain dependence on the size parameter A^. 

For the orthogonal polynomials Pn{x) satisfying the recursion relation 



there exits a powerful theorem by Van Doorn on the lower bound on the lowest zero 
which has been generalised by Ismail and Li ^ to include an upper bound to the largest 
zero. Simply stated Ismail and Li's result is the following: if x^j}^ and x^j^'' are the smallest 
and largest zeros respectively of Pk{x) with A; > 1 then they are bounded by the interval 
{A, B) where 



Dn{x) = {an - x)Dn-l{x) - (3l_]^Dn-2{x) 



(11) 



Pn{x) = {X- an)Pn-l{x) - 0^_iPn-2 



(12) 



A = min{/„ : 1 < n < k} 



B = max{/^ : 1 < n < k} 



(13) 



and the bound sequence is given by 




(14) 
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In the above expression for f^, {an}i ^ is a chain sequence. That is, there exists a parameter 
sequence {gn}i for which we have the factorisation of a„ 

an = 9n+iii- 9n), l<n<k (15) 

where < gi < 1 and < < 1 for 1 < n < k. Van Doom's theorem is for strict 
equahty, whereby maximising the lower bound with respect to the parameter sequence gives 
the lowest eigenvalue exactly: 



(1) 

xi = max 
{9} 



min{/^ : 1 < n < k} 



(16) 



for k > 1. Only a technical difference occurs with this theorem in that gi = and gk = 1- 

For finite the termination of the Lanczos recursion occurs at some nmax when the 
sector of Hilbert space, determined by the trial state, has been exhausted. In general the 
basis of states grows faster than any linear enumeration with A^. Taking the minimum of 
to occur at some n = zN [n < nmax), the asymptotic forms for ^^(A^) and PniN) give 
the leading order behaviour of the lowest energy level as 



x^^ = max 



'-max 



aiz)- ^(3{z) + 0{l/N) 



(17) 



{9} 

Since (5n{N) > 0, to maximise the RHS of the above expression with respect to the 
chain sequence we may choose the maximal constant chain sequence |Q 

a = ^ ^. (18) 

4cos^ 



We can now obtain the ground state energy density, Sq, in the bulk limit: 

£0= lim xi^) . (19) 

In the N —>■ 00 limit, we also have rimax 00 giving a — » 1/4, and it is straightforward 
to establish that in general z is obtained by finding the greatest lower bound with respect 
to z. The ground state energy density of the general extensive many-body problem in the 
bulk limit is therefore 

£o = inf [«(^)-2/5(^)]. (20) 
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IV. EXAMPLES 



As a first illustration of this exact analytic diagonalisation of extensive systems we 
consider the purely mathematical model defined by the tri-diagonal form 



where r = a/b. The minimum occurs at z^^^ = ln(2/r)/6(r — 1). It is a straightforward 
matter to numerically diagonalise this system for increasing to demonstrate the conver- 
gence of the numerical results to the analytic expression. The results of this exercise are 
shown in Figure 1 (for a typical case a = 3 and 6 = 2) where we plot the error, defined 
as the difference from the exact value, of numerical diagonalisations for increasing A^. The 
approach to the analytic solution is clear as ^ cxd. 

Secondly we compare the the analytic diagonalisation with numerical results using the 
truncated Lanczos coefficients given by the plaquette expansion. We do this here for the 
case of the anti-ferromagnetic Heisenberg chain, for which the expansions have been derived 
to high order using the classical Neel state as the trial state. This state is an poor choice 
for the isotropic Heisenberg model but is sufficiently simple to allow the computation of 
moments up to < H^^ >c- Employing the traditional analysis of the plaquette expansion 
the Tn matrices are diagonalised numerically for increasing chain size f^. In Figures 2 
and 3 we demonstrate the approach of the numerical diagonalisation for the ground state 
to the analytic expression as the N oo limit is approached. We show the two typical 
cases which can occur due to the truncation of the expansions for a{z) and l3{z). Figure 2 
corresponds to an order of the plaquette expansion for which a{z) — 2(3{z) has no minimum 
- a point of inflection develops, reflecting the fact that the expansion naturally breaks down 




(21) 



The lowest eigenvalue in the N oo limit is 
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at some large enough value of z. Figure 3 shows a case where a minimum develops - the 
numerical values for the lowest eigenvalue in the large volume limit match onto the value of 
at the minimum. 

We must emphasis that the plaquette expansion is a series expansion of the exact 
Lanczos coefficients about 2 = and truncated, and is therefore only reliable for small z. 
The error involved grows rapidly with z which can occur with either sign, and thus the 
large z behaviour of the truncated approximations bear no relation to the exact behaviour. 
Because of this we don't expect that at every order at which the truncation is made, that a 
minima would arise. 

V. CONCLUSIONS 

We have found an expression for the ground state energy density for the extensive 
many-body problem which is completely general and, since the Hamiltonian was diagonalised 
exactly, is non-perturbative. Given the Hamiltonian in tri-diagonal form the expression can 
be used immediately. So far, the exact analytic transformation of a system to tri-diagonal 
form has not been achieved for any examples of solvable systems. This state of affairs may 
change as the theorem proved here gives impetus to efforts in this direction. However, an 
immediate approximate tri-diagonalisation of the general problem does exist in the plaquette 
expansion and studies of the usefulness of this method with better trial states also merits 
further work. The analytic nature of this expansion uncovered the existence of the scaled 
Lanczos iteration parameter z = n/N and the asymptotic forms a{z) and (3{z) for the exact 
problem. 

An interesting extension of this work is to consider a similar analysis on excited states. 
An indication that the solution for the mass-gap of the general extensive many-body system 
may in fact be possible is the fact that the mass-gap of the first order plaquette expansion 
has already been solved analytically Such a generalisation is presently under active 

investigation. 
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FIGURES 

FIG. 1. Comparison of the analytic solution with numerical diagonalisation for the simple case, 
a„(iV) = (1 + an/N^)^ and /3n(iV) = (1 + bn/N^)^ - 1, with a = 3 and 6 = 2. 

FIG. 2. Plaquette expansion of the ID AFH model at an order where a point of 

inflection develops due to the breakdown of the expansion. The numerical diagonalisation for 
increasing N converges to the function a{z) — 2(3{z). 

FIG. 3. Plaquette expansion of the ID AFH model at an order (l/AT^) where a minimum 
develops in a{z) — 2f3{z). The numerical diagonalisation data &i N = b x 10"^ clearly display the 
convergence to the value of min[Q;(2;) — 2(3{z)\. 
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